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I. INTRODUCTION 

Two-dimensional (2D) Y-electron systems (with a 
small finite N) in strong magnetic fields (B) have been 
the focus of extensive theoretical investigations in the last 
twenty years, 1 ! 2 . 3 . 4 . 5 . 6 . 7 . 8 . 9 . 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 The princi- 
pal motivations for these research activities are: (I) The 
early realization 12 that certain special states of few elec- 
tron systems are relevant* 1 ^ through appropriate analo- 
gies to the physics of the fractional quantum Hall ef- 
fect (FQHE), observed in the infinite 2D electron gas; 
(11) The unavoidable necessity, due to computer limi- 
tations, to test proposed model wave functions for the 
FQHE through numerical calculations for finite-size sys- 
tems; (III) The recent progress in nanofabrication tech- 
niques at semiconductor interfaces that has allowed ex- 
periments on 2D quantum dots (QD's), with refined con- 
trol of their size, shape, and number of electronsi2i22i2i 
(down to a few electrons). 

The physics of such systems (i.e., QD's in high B), 
is most often described with the use of composite- 
fermion 4 / Jastrow-Laughlin 1 (CF/JL) analytic trial wave 
functions in the complex plane. However, it is well known 
that the thematic framework of the CF/JL approach 
is built on the so-called Jastrow correlations associated 
with a particular short-range interparticle repulsion^ In 
a recent paper^ using as a thematic basis the picture 
of collectively rotating electron (or Wigner) molecules 
(REM's), we have derived a different class of analytic 
and parameter-free trial wave functions. The promising 
property of these REM wave functions is that, unlike the 
CF/JL ones, they capture the all-important correlations 
arising from the long-range character of the Coulomb 
force. 

In this paper, we present an in-depth assessment of 
the CF/JL and REM trial wave functions regarding their 
ability to approximate the exact wave functions in the 
case of QD's (this case is often referred to as the "disk 
geometry" in the FQHE literature) . First systematic ex- 
act diagonalization (EXD) results are reported here for 
the lowest rotational band of N = 6 electrons in strong 



magnetic fields in the range of high angular momenta 
70 < L < 140 (covering the corresponding range of 
fractional filling factors^ 3 . 1/5 > v > 1/9). A detailed 
comparison (addressing five properties, i.e., prediction of 
magic angular momenta, radial electron densities, occu- 
pation number distributions, overlaps and total energies, 
and exponents of current-voltage power law) shows that 
the REM wave functions yield a superior description to 
that obtained through the composite-fermion/ Jastrow- 
Laughlin ones. 

The plan of this paper is as follows: Section II presents 
an outline of the REM theory, while section III focuses on 
a brief review of the composite-fermion approach. Exact- 
diagonalization results and comparisons with the CF/JL 
and REM wave functions are presented in section IV. 
Finally, our results are summarized in section V. 



II. OUTLINE OF REM THEORY 

In the last eight years, and in particular since 
1999 [when it was demonstrated 2 ^ that Wigner crys- 
tallization is related to symmetry breaking at the 
unrestricted Hartree-Fock (UHF) mean-field level], 
the number of publications 8-17 ' 24-40 addressing the 
formation and properties of Wigner (or electron) 
molecules in 2D QD's and quantum dot molecules 
has grown steadily. A consensus has been reached 
that rotating electron molecules are formed both 

m z „ r 12. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34.35. 36. 37. 38. 39. 40 an( j 

high 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 magnetic fields. 

At B = 0, formation of REM's in QD's is analogous 
to Wigner crystallization in infinite 2D media, i.e., when 
the strength of the interelectron repulsion relative to the 
zero-point kinetic energy (Rw) exceeds a certain criti- 
cal value, electrons spontaneously crystallize around sites 
forming geometric molecular structures. At high mag- 
netic fields, the formation of Wigner molecules may be 
thought of as involving a two-step crystallization process: 
(I) the localization of electrons results from the shrinkage 
of the orbitals due to the increasing strength of the mag- 
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netic field; (II) then, even a weak interelectron Coulomb 
repulsion is able to arrange the localized electrons accord- 
ing to geometric molecular structures (thus this process is 
independent of the value of Rw)- It has been foun d 8 ' 10 i 12 
that the molecular structures at high B coincide with the 
equilibrium configurations at B = of N classical point 
charges 41ti£ 

Due to the finite number, N, of electrons, however, 
there are two crucial differences between the REM and 
the bulk Wigner crystal. Namely, (I) the crystalline 
structure is that of the equilibrium 2D configuration of 
N classical point charges, and thus consists of nested 
polygonal rings;— (II) the Wigner molecules rotate as a 
whole (collective rotations) in analogy with the case of 
3D natural molecules. 

A most striking observation concerning the REM's 
is that their formation and properties have been es- 
tablished with the help of traditional ab initio many- 
body methods, i.e., exact diagonalization^i2iii*iSiS42L2£ 
quantum Monte Carlo 26 i 29 i 33 i 39 (QMC), and the system- 
atic controlled hierarchj&iSiiSiiiiSiMiSiiii of approxima- 
tions involving the UHF and subsequent post-Hartree- 
Fock methods. This contrasts with the case of the CF/JL 
wave functions, which were inspired through "intuition- 
based guesswork" . 

In spite of its firm foundation in many-body theory, 
however, the REM picture has not, until recently, suc- 
cessfully competed with the CF/JL picture; indeed many 
research p a p e rsi^^S4LiS42iS2iSi and booksi& describe 
the physics of QD's in high magnetic fields following 
exclusively notions based on CF/JL functions, as ex- 
pounded in 1983 (see Ref. Q) and developed in detail in 
1995 in Ref. and Ref. 0. We believe that one of the 
main obstacles for more frequent use of the REM picture 
has been the lack of analytic correlated wave functions 
associated with this picture. This situation, however, has 
changed with the recent explicit derivation of such REM 
wave functions— 

The approach used in Ref. for constructing the 
REM functions in high B consists of two-steps: First the 
breaking of the rotational symmetry at the level of the 
single-determinantal unrestricted Hartree-Fock approxi- 
mation yields states representing electron molecules (or 
finite crystallites, also referred to as Wigner molecules, 
see Ref.l24l and Ref. 0). Subsequently the rotation of 
the electron molecule is described through restoration 
of the circular symmetry via post Hartree-Fock meth- 
ods, and in particular Projection Techniques— Natu- 
rally, the restoration of symmetry goes beyond the single 
determinantal mean-field description and yields multi- 
determinantal wave functions. For QD's, we have shown 
that the method of symmetry restoration is applicable to 
both the zero^^ifii and high— magnetic-field cases. 

In the zero and low-field cases, the broken symmetry 
UHF orbitals need to be determined numerically, and, 
in addition, the restoration of the total-spin symmetry 
needs to be considered for unpolarized and partially po- 
larized cases. The formalism and mathematical details 



of this procedure at B = have been elaborated in Ref. 
l34l (see also Ref. |S3 and Ref. l54|) for the restoration of 
the total spin in the case of quantum dot molecules). 

In the case of high magnetic fields, one can specifi- 
cally consider the limit when the confining potential can 
be neglected compared to the confinement induced by 
the magnetic field. Then, assuming a symmetric gauge, 
the UHF orbitals can be representee^*^ by displaced 
Gaussian analytic functions, centered at different posi- 
tions Zj = Xj + iYj according to the equilibrium config- 
uration of N classical point charges^i*^ arranged at the 
vertices of nested regular polygons (each Gaussian rep- 
resenting a localized electron) . Such displaced Gaussians 
are written as (here and in the following 1 = >/— 1) 

u(z , Z j ) = (l/V^) 

x expl-lz-Ztffflexpl-iixYj-yXj)], (1) 

where the phase factor is due to the gauge invariance. 
z = x + ly (see Ref. IHfih . and all lengths are in dimen- 
sionless units of IbV% with the magnetic length being 
Ib = y/hc/eB. 

In Ref. ITHl we used these analytic orbitals to first con- 
struct the broken symmetry UHF determinant, , J'^, HF , 
and then proceeded to derive analytic expressions for 
the many-body REM wave functions by applying onto 
^P^ HF an appropriate projection operator— Ol that re- 
stores the circular symmetry and generates correlated 57 
wave functions with good total angular momentum L. 
These REM wave functions can be easily written downi& 
in second-quantized form for any classical polygonal 
ring arrangement (ni,ri2, ...) by following certain simple 
rules for determining the coefficients of the determinants 
D(li,l 2 , —,In) = det[z^, z l 2, ■ ■ •, z$], where the l/s de- 
note the angular momenta of the individual electrons. 
Since we will focus here on the case of N = 6 and N = 3 
electrons, we list for completeness the REM functions as- 
sociated with the (0, N) and (1, N—l) ring arrangements, 
respectively [here (0, N) denotes a regular polygon with 
N vertices, such as an equilateral triangle or a regular 
hexagon, and (1, N — 1) is a regular polygon with N — 1 
vertices and one occupied site in its center], 
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L = L + (N - \)m, m = 0,1,2,3,. 
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where Lq = N(N — l)/2 is the minimum allowed total 
angular momentum for (polarized) electrons in high 
magnetic fields. 

Notice that the REM wave functions [Eq. © and Eq, 
(@J] vanish identically for values of the total angular mo- 
menta outside the specific values given by Eq. @ and 
Eq. (JS}, respectively. 



III. OUTLINE OF COMPOSITE-FERMION 
THEORY 

According to the CF picture^ the many body wave 
functions in high magnetic fields that describe TV- 
electrons in the disc geometry (case of 2D QD's) are given 
by the expression, 
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where z = x+iy and is the Slater determinant of N 

non-interacting electrons of total angular momentum L*; 
it is constructed according to the Independent Particle 
Model (IPM) from the Darwin- FockS orbitals ip P! i(z), 
where p and I are the number of nodes and the angular 
momentum, respectively [for the values of p and I in the 
nth Landau level in high B, see the paragraph following 
equation J7J below]. 

The Jastrow factor in front of , J^ M is introduced to 
represent the effect of the interelectron Coulombic inter- 
action. In the CF literature, this assumption is often 
described by saying that "the Jastrow factor binds 2m 
vortices to each electron of <I'^ M to convert it into a 
composite fermion" . 

The single-particle electronic orbitals in the Slater de- 
terminant , J^ M are not restricted to the lowest Landau 
level (LLL). As a result, it is necessary to apply a pro- 
jection operator Plll to guarantee that the CF wave 
function lies in the LLL, as appropriate for B — > oo. 

Since the CF wave function is an homogeneous poly- 
nomial in the electronic positions Zj's, its angular mo- 
mentum L is related to the non-interacting total angular 
momentum L* as follows, 



L* + mN(N - 1) = L* + 2mL . 



(7) 



There is no reason to a priori restrict the Slater de- 
terminants ^Y' M to a certain form, but according to 
Ref. H such a restriction is absolutely necessary in or- 
der to derive systematic results. Thus following Ref. 
6, henceforth, we will restrict the non-interacting L* to 
the range — Lq < L* < Lq, and we will assume that 
the Slater determinants ^5> l £* M are the so-called com- 
pact ones. Let N n denote the number of electrons in 
the nth Landau Level (LL) with J2n=o N n = N; t is 
the index of the highest occupied LL and all the lower 
LL's with n < t are assumed to be occupied. The 
compact determinants are defined as those in which the 
N n electrons occupy contiguously the single-particle or- 
bitals (of each nth LL) with the lowest angular momenta, 
/ = -n, -n + 1, — n + N n — 1 [p + - l)/2 = n]. 
The compact Slater determinants are usually denoted as 
[No, Ni, N t ], and the corresponding total angular mo- 
menta are given by L* = (1/2) Y? s =o N *( N s - 2s - 1). 

Most important for our present study is the fact that 
the Jastrow-Laughlin wave functions with angular mo- 
mentum L — (2m + 1)Lq [corresponding to fractional 
filling factors v = L /L = l/(2m + 1)], 

$t h (N)= J] fe-%) 2ro+1 e X p(-^^/2), 

l<i<j<N \ fc=l / 

(8) 

are a special case of the CF functions for L* — Lq, i.e., 

$i L (A0 = $£ F (iV;L* = Lo), L=(2m+l)L . (9) 

Note that for L* = Lq, all the non- interacting electrons 
occupy contiguous states in the LLL (n = 0) with I = 
Q,1,..,AT-1. 

The CF/JL wave functions [equations © and iJSJ] are 
represented by compact, one-line mathematical expres- 
sions, which however are not the most convenient for 
carrying out numerical calculations. Numerical studies of 
the CF/JL functions usually employ sophisticated Monte 
Carlo computational techniques. The REM wave func- 
tions, on the other hand, are by construction expressed in 
second-quantized (superposition of Slater determinants) 
form, precisely like the wave functions from exact diag- 
onalization, a fact that greatly simplifies the numerical 
work. In the numerical calculations involving JL wave 
functions in this paper, we have circumvented the need 
to use Monte Carlo techniques, since we were able to de- 
termine the Slater decomposition 59 of the JL states with 
the help of the symbolic language MATHEMATICAL 

We stress again that, unlike the REM functions, the 
CF / JL wave functions have not been derived microscop- 
ically, i.e., from the many-body Schrodinger equation 
with interelectron Coulombic repulsions. Attempts have 
been made to justify them a posteriori by pointing out 
that their overlaps with exact wave functions are close to 
unity or that their energies are close to the exact ener- 
gies. However, we will show below that this agreement 
is limited to rather narrow ranges of filling factors be- 
tween 1 > v > 1/3 or to small electron numbers N; as 
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soon as one extends the comparisons to a broader range 
of u's for N > 6, as well as to other quantities like elec- 
tron densitiess and occupation number distributions, this 
agreement markedly deteriorates. 



IV. EXACT DIAGONALIZATION RESULTS 
AND COMPARISONS 

In the case of high magnetic fields, the Hilbert space 
for exact-diagonalization calculations can be restricted 
to the LLL and many such calculations have been 
rcported 2 i 3 i 5 i 9 i 10 i 18 . 44 i 46 i 49 i 61 i 62 i 63 i 64 in the past twenty 
years. However, for N > 5, such EXD studies have 
been restricted to angular momenta corresponding to the 
rather narrow range of fillings factors 1 > v > 1/3. 

In this paper, we have performed systematic EXD cal- 
culations in the LLL for N = 6 electrons covering the 
much broader range of fillings factors 1 > v > 1/9; such 
a range corresponds to angular momenta 15 < L < 140 
(note that for v = 1/3 one has L = 45). Of crucial im- 
portance for extending the calculations to such large L's 
has been our use of Tsiper'sS 5 - analytic formula for cal- 
culating the two-body matrix elements of the Coulomb 
intcrclcctron repulsion; this formula expresses the ma- 
trix elements as finite sums of positive terms. Earlier 
analytic formulas 3 suffered from large cancellation errors 
due to summations over alternating positive and negative 
terms. At the same time, Tsiper's formula is computa- 
tionally faster compared to the slowly-convergent series 
of Ref . EH 

For the solution of the large scale, but sparse, Coulomb 
eigenvalue problem, we have used the ARPACK com- 
puter codei£& For a given L, the Hilbert space is built 
out of Slater determinants, 

N 

L)(Z 1 ,Z 2 ,..,^)e X p(-^ i z*/2), ( 10 ) 

i=l 



A. Predictions of magic angular momenta 

Foe iV = 6, Figs. 1-4 display (in four installments) the 
total interaction energy from EXD as a function of the 
total angular momentum L in the range 19 < L < 140. 
(The total kinetic energy, being a constant, can be dis- 
regarded.) One can immediately observe the appearance 
of downward cusps, implying states of enhanced stability, 
at certain "magic angular momenta" . 

For the CF theory, the magic angular momenta can be 
determined by Eq. J7J, if one knows the non- interacting 
L*'s; the CF magic L's in any interval l/(2m — 1) > 
v > l/(2m + 1) [15(2m - 1) < L < 15(2m + 1)], m = 
1,2,3,4, can be found by adding 2mLo = 30m units 
of angular momentum to each of the L*'s. To obtain the 
non-interacting L*'s, one needs first to construd**^ the 
compact Slater determinants. The compact determinants 
and the corresponding non-interacting L*'s are listed in 
Tabled 

There are nine different values of L*'s, and thus the 
CF theory for N = 6 predicts that there are always nine 
magic numbers in any interval 15(2m— 1) < L < 15(2m+ 
1) between two consecutive JL angular momenta 15(2m— 
1) and 15(2m + 1), m = 1,2,3,... (henceforth we will 
denote this interval as T m ). For example, using Tabled 
and Eq. Ql, the CF magic numbers in the interval 15 < 
L < 45 (m = 1) are found to be the following nine^i 

15, 21, 25, 27, 30, 33, 35, 39, 45. (12) 

On the other hand, in the interval 105 < L < 135 (m = 
4), the CF theory predicts the following set of nine magic 
numbers, 

105, 111, 115, 117, 120, 123, 125, 129, 135. (13) 

An inspection of the total-energy- vs. -L plots in Figs. 
1-4 reveals that the CF prediction badly misses the ac- 
tual magic angular momenta specified by the EXD cal- 
culations as those associated with the downward cusps. 



with 

N 

h <l 2 < ... <l N , ^2h=L, (li) 
fc=l 

and its dimensions are controlled by the maximum al- 
lowed single-particle angular momentum l max , such that 
h < 'ma, 1 < k < N. We have used Z max = Z^ L ax + 5 = 
10(to+ 1) (see Ref. |5!|for the definition of /,^ x ) for each 
group of angular momenta L corresponding to the range 
l/(2m + 1) < v < l/(2m - 1), m = 1,2,3,4. For ex- 
ample, for L = 105, i max = 40 and the dimension of the 
Hilbert space is 56115; for L = 135, Z max = 50 and the 
size of the Hilbert space is 187597. By varying Z max , we 
have checked that this choice produces well converged 
numerical results. 



TABLE I: Compact non-interacting Slater determinants and 
associated angular momenta L* for N — 6 electrons accord- 
ing to the CF presciption. Both L* — —3 and L* = 3 are 
associated with two compact states each, the one with lowest 
energy being the preferred one. 



Compact state 


L* 


[1,1,1,1,1,1] 


-15 


[2,1,1,1,1] 


-9 


[2,2,1,1] 


-5 


[3,1,1,1] 


-3 


[2,2,2] 


-3 


[3,2,1] 





[4,1,1] 


3 


[3,3] 


3 


[4,2] 


5 


[5,1] 


9 


[6] 


15 
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FIG. 1: Total interaction energy from exact-diagonalization 
calculations as a function of the total angular momentum 
(10 < L < 50) for N = 6 electrons in high magnetic field. The 
upwards pointing arrows indicate the magic angular momenta 
corresponding to the classically most stable (1,5) polygonal 
ring arrangement of the Wigner molecule. The short down- 
wards pointing arrows indicate successful predictions of the 
composite-fermion model. The long downward arrow indi- 
cates a magic angular momentum not predicted by the CF 
model. Energies in units e 2 /kIb, where k is the dielectric 
constant. 



Indeed it is immediately apparent that the number of 
downward cusps in any interval X m is always different 
from 9. Indeed, there are 10 cusps in I\ (including that 
at L = 15, not shown in Fig.QJ, 10 in^ (see Fig.|2J), 7 in 
1 3 (see Fig. EJ|, and 7 in X 4 (see Fig.gJ). In detail, the CF 
theory fails in the following two aspects: (I) There are 
exact magic numbers that are consistently missing from 
the CF prediction in every interval; with the exception of 
the lowest L = 20, these exact magic numbers (marked 
by a long downward arrow in the figures) are given by 
L = 10(3to- 1) and L = 10(3™ + 1), to = 1,2,3,4,...; 
(II) There are CF magic numbers that do not correspond 
to downward cusps in the EXD calculations (marked by 
medium-size downward arrows in the figures). This hap- 
pens because cusps associated with L's whose difference 
from Lq is divisible by 6 (but not simultaneously by 5) 
progressively weaken and completely disappear in the in- 
tervals T m with m > 3; only cusps with the difference 
L — Lq divisible by 5 survive. On the other hand, the 
CF model predicts the appearance of four magic num- 
bers with L — Lq divisible solely by 6 in every interval 
T m , at L = 30m T- 9 and 30to =f 3, m = 1, 2, 3, ... The 
overall extent of the inadequacy of the CF model can 
be appreciated better by the fact that there are six false 
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FIG. 2: Total interaction energy from exact-diagonalization 
calculations as a function of the total angular momentum 
(40 < L < 80) for N — 6 electrons in high magnetic field. The 
upwards pointing arrows indicate the magic angular momenta 
corresponding to the classically most stable (1,5) polygo- 
nal ring arrangement of the Wigner molecule. The short 
downwards pointing arrows indicate successful predictions of 
the composite-fermion model. The medium-size downwards 
pointing arrow indicates a prediction of the CF model that 
fails to materialize as a magic angular momentum. The long 
downward arrows indicate magic angular momenta not pre- 
dicted by the CF model. Energies in units of e 2 /kIb, where 
k is the dielectric constant. 



predictions (long and medium-size downward arrows) in 
every interval I m with to > 3, compared to only five cor- 
rect ones (small downward arrows, see Fig. |3 and Fig. 
EJ). 

In contrast to the CF model, the magic angular mo- 
menta in the REM theory are associated with the polygo- 
nal ring configurations of N classical point charges. This 
is due to the fact that the enhanced stability of the down- 
ward cusps results from the coherent collective rotation of 
the regular-polygon REM structures. Due to symmetry 
requirements, such collective rotation can take place only 
at magic-angular-momenta values. The in-between angu- 
lar momenta require the excitation of additional degrees 
of freedom (like the center of mass and/or vibrational 
modes) , which raises the total energy with respect to the 
values associated with the magic angular momenta. 

For N = 6, the lowest in energy ring configuration 
is the (1,5), while there exists a (0,6) isomer^i*^ with 
higher energy. As a result, our EXD calculations (as well 
as earlier onesiii*ii for lower angular momenta L < 70) 
have found that there exist two sequences of magic an- 
gular momenta, a primary one (S p ) with L = 15 + 5to 
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FIG. 3: Total interaction energy from exact-diagonalization 
calculations as a function of the total angular momentum 
(70 < L < 110) for N = 6 electrons in high magnetic 
field. The upwards pointing arrows indicate the magic an- 
gular momenta corresponding to the classically most stable 
(1,5) polygonal ring arrangement of the Wigner molecule. 
The short downwards pointing arrows indicate successful pre- 
dictions of the composite-fermion model. The medium-size 
downwards pointing arrows indicate predictions of the CF 
model that fail to materialize as magic angular momenta. The 
long downward arrows indicate magic angular momenta not 
predicted by the CF model. Energies in units of e 2 /kIb, where 
k is the dielectric constant. 



FIG. 4: Total interaction energy from exact-diagonalization 
calculations as a function of the total angular momentum 
(100 < L < 140) for N = 6 electrons in high magnetic 
field. The upwards pointing arrows indicate the magic an- 
gular momenta corresponding to the classically most stable 
(1,5) polygonal ring arrangement of the Wigner molecule. 
The short downwards pointing arrows indicate successful pre- 
dictions of the composite-fermion model. The medium-size 
downwards pointing arrows indicate predictions of the CF 
model that fail to materialize as magic angular momenta. The 
long downward arrows indicate magic angular momenta not 
predicted by the CF model. Energies in units of e 2 /kIb, where 
k is the dielectric constant. 



[see Eq. JJJJ], associated with the most stable (1,5) classi- 
cal molecular configuration, and a secondary one (S s ) 
with L = 15 + 6m [see Eq. ||SJ|], associated with the 
metastable (0, 6) ring arran gemen t. Furthermore, our 
calculations (see also Refs. Illll.Sh show that the sec- 
ondary sequence S s contributes only in a narrow range of 
the lowest angular momenta; in the region of higher an- 
gular momenta, the primary sequence S p is the only one 
that survives and the magic numbers exhibit a period of 
five units of angular momentum. It is interesting to note 
that the initial competition between the primary and sec- 
ondary sequences, and the subsequent prevalence of the 
primary one, has been seen in other sizes as well^i i.e., 
N = 5,7, 8. Furthermore, this competition is reflected 
in the field-induced molecular phase transitions associ- 
ated with broken symmetry UHF solutions in a parabolic 
QD. Indeed, Ref. lljj demonstrated recently that, as a 
function of increasing B, the UHF solutions for TV = 6 
first depict the transformation of the maximum-density- 
drople^ into the (0,6) molecular configuration; then (at 
higher B) the (1,5) configuration replaces the (0,6) struc- 
ture as the one having the lower HF energy^ 



The extensive comparisons in this subsection lead in- 
evitably to the conclusion that the CF model cannot ex- 
plain the systematic trends exhibited by the magic angu- 
lar momenta in 2D QD's in high magnetic fields. These 
trends, however, were shown to be a natural consequence 
of the formation of REM's and their metastable isomers. 



B. Radial electron densities 

We turn now our attention to a comparison of the ra- 
dial electron densities (ED's). Fig. displays the corre- 
sponding ED's from EXD, REM, and CF/JL wave func- 
tions at three representative total angular momenta, i.e., 
L = 75 (v = 1/5), 105 (1/7), and 135 (1/9). 

An inspection of Fig. \E\ immediately reveals that (I) 
The EXD radial ED's (solid lines) exhibit a prominent 
oscillation corresponding to the (1,5) molecular struc- 
ture (averaged over the azimuthal angles). Indeed the 
integral of the exact ED's from the origin to the mini- 
mum point between the two humps is practically equal 
to unity; (II) There is very good agreement between 
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FIG. 5: Radial electron densities for TV = 6 electrons in 
high magnetic field. Solid line: densities from exact diago- 
nalization. Dashed line: densities from REM wave functions. 
Dotted line: densities from Jastrow-Laughlin wave functions. 



FIG. 6: Distribution of occupation numbers as a function 
of single-particle angular momentum I for N = 6 electrons 
in high magnetic field. Solid circles: occupation numbers 
from exact diagonalization. Open circles: occupation num- 
bers from REM wave functions. Crosses: occupation numbers 
from Jastrow-Laughlin wave functions. 



the REM (dashed lines) and exact ED's; this agreement 
improves with higher angular momentum; (III) The JL 
ED's (dotted lines ) miss the oscillation of the exact ED 
in all three cases in a substantial way. 

The inability of the radial ED's calculated with the JL 
functions to capture the oscillations exhibited by the ex- 
act ones was also seen recently for the v = 1 /3 case and 
for all electron numbers iV = 6, 7, 8, 9, 10, 11, 12 in Ref.[H 
(see in particular Fig. 1 therein). We further note that 
the oscillations of the exact ED's in that figure corre- 
spond fully to the classical molecular ring arrangements 
listed in Ref. |U e.g., to (1,7) for N = 8 and to (3,9) 
for N — 12, in agreement with our rotating-electron- 
molecule interpretation. 



C. Distribution of occupation numbers 

In this subsection, we address the behavior of the 
occupation-number distribution n(l) = ($|a|a;|$) as 
a function of the single-particle angular momentum I, 
where the creation and annihilation operators refer to 
the single-electron states ipo t i(z) in the LLL. For N = 6, 
Fig. displays the n(/)'s from all three families of wave 
functions, i.e., EXD (solid circles), REM (open circles), 
and JL (crosses), and for the three representative angular 
momenta L = 75 [y = 1/5), 105 (1/7), and 135 (1/9). 

Again, an inspection of Fig.^immediately reveals that 
(I) The EXD occupation numbers exhibit a prominent os- 
cillation corresponding to the (1,5) molecular structure. 
Indeed the sum of the exact n(l)'s from I = to the min- 
imum point between the two humps is practically equal 
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TABLE II: Case of N = 3 electrons in high mag- 
netic fileds. Overlaps, ($ £ |$f x )/«$z, <$! x |$! x » 1/2 , of 
REM's (<&'s) and JL functions ($'s) with the corresponding 
exact eigenstates ($ EX 's) for various values of the angular mo- 
menta L (y are the corresponding fractional filling factors). 
Recall that the angular momenta for the JL functions are 
L JL = N(N ~ l)(2m + l)/2, with m = 0, 1, 2, 3, ... The JL 
overlaps are from Ref. [lj 



L{v) JL REM 

9(1/3) 0.99946 0.98347 

15(1/5) 0.99468 0.99473 

21(1/7) 0.99476 0.99674 

27(1/9) 0.99573 0.99758 

33(1/11) 0.99652 0.99807 

39(1/13) 0.99708 0.99839 



to unity; (II) There is very good agreement between the 
REM and exact occupation numbers; this agreement im- 
proves with higher angular momentum; (III) For all three 
cases, the JL occupation numbers exhibit a systemati- 
cally different trend and they are not able to capture the 
oscillatory behavior of the EXD occupation numbers. 

We further note that a substantial discrepancy between 
JL and EXD occupation numbers was also noted in Ref. 
for the case of N = 7 electrons and v = 1/3 (L = 63). 

The systematic deviations between the JL and EXD 
ED's and occupation numbers inevitably points to the 
conclusion that these two families of wave functions rep- 
resent very different many-body physical problems. In- 
deed, the JL functions have been found^ to be exact so- 
lutions for a special class of short-range two-body forces, 
while the EXD functions faithfully reflect the long-range 
character of the Coulombic interelectron repulsion. On 
the other hand, as discussed in Ref. 15, the REM wave 
functions, derived through a traditional many-body ap- 
proach, are able to capture the correlations arising from 
the long-range character of the Coulomb force; the os- 
cillatory behavior of the EXD and REM ED's and oc- 
cupation numbers (associated with formation of Wigner 
molecules) constitutes a prominent and unmistaken sig- 
nature of such Coulombic correlations. 



TABLE III: Overlaps of JL and REM wave functions with the 
exact ones for N = 6 electrons and various angular momenta 
L (y are the corresponding fractional filling factors). 





JL 


REM 


75(1/5) 


0.837 


0.817 


105(1/7) 


0.710 


0.850 


135(1/9) 


0.665 


0.860 



D. Comparison of overlaps and total energies 

We turn now our attention to the overlaps of the REM 
and JL wave functions with those obtained through ex- 
act diagonalization. We start by listing in Table [H] the 
overlaps for the simpler case of N — 3 electrons in high 
magnetic fields. One sees immediately that these over- 
laps are all very close to unity (> 0.99) for both the REM 
and JL cases and for even rather high angular momenta 
[e.g., L = 39 {v = 1/13)]. 

Ever since they were calculated by Laughlin in his orig- 
inal paper the JL overlaps for N = 3 electrons have ex- 
ercised a great influence in the literature of the fractional 
quantum Hall effect (FQHE) . Indeed, in a rather sweep- 
ing generalization to any N and L (note that Ref. ^ has 
indeed found that the JL overlaps for v = 1/3 remain 
very close to unity for all cases with 5 < N < 12), the 
close-to-unity values of the JL overlaps have been pre- 
sumed to provide "proof" that the CF/JL functions ap- 
proximate very well the corresponding exact many-body 
wave functions; as we have already shown earlier, this 
presumption is highly questionable. 

We have calculated the overlaps for N = 6 elec- 
trons and for the three representative higher angular- 
momentum values L — 75 (v — 1/5), 105 (1/7), and 
135 (1/9); the results are listed in Table ITTT1 for both the 
REM and JL wave functions. A most remarkable feature 
of the results in Table inTl is that the extraordinary, higher 
than 0.99 values (familiar from Laughlin' s papesi) are 
totally absent. Instead, the JL overlaps rapidly deteri- 
orate for higher L's (lower v's), and for v = 1/9 they 
have attained values below 0.67. In contrast, the REM 
overlaps remain above 0.80 and slowly approach unity as 
L increases. 

From our results for v < 1/5 and the results of Ref. 
1621 for v = 1/3, it is apparent that the overlaps alone 
are not a reliable index for assessing the agreement or 
disagreement between trial and exact wave functions. For 
example, for N = 6 and L = 75 (y = 1/5), Table llTTl 
shows that the JL and REM overlaps are close to each 
other (0.837 vs. 0.817). However, as the earlier analyses 
based on the electron densities and occupation numbers 
show, the JL wave function is not a good approximation 
to the exact one; in contrast, the REM wave function 
offers a much better description. 

In addition to the overlaps, earlier studies (see, e.g., 
Ref. El have also relied on the total energies for assess- 
ing the agreement, or not, between CF and exact wave 
functions. We thus list in Table ITVl the total energies for 
N = 6 and for the three representative higher angular- 
momentum values L = 75 {y — 1/5), 105 (1/7), and 
135 (1/9). It is seen that both the JL and REM total 
energies exhibit very small relative errors compared to 
the corresponding EXD ones in all three instances, a fact 
that indicates that, by themselves, the total energies^ 
are an even less reliable index compared to the overlaps. 
In particular, note that for N = 6 and L = 135, the JL 
and exact total energies differ only in the third decimal 
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TABLE IV: Total interaction energies of JL, REM, and exact- 
diagonalization wave functions for TV = 6 electrons and var- 
ious angular momenta L (y are the corresponding fractional 
filling factors). The percentages within parentheses indicate 
relative errors. Recall that the angular momenta for the JL 



functions are L JL = N(N - l)(2m + l)/2, m = 0,1,2,3,... 
Energies in units of e 2 /kIb, where k is the dielectric constant. 

L{v) JL REM EXACT 

75(1/5) 2.2093 (0.32%) 2.2207 (0.85%) 2.2018 

85(3/17) 2.0785 (0.65%) 2.0651 

95(3/19) 1.9614 (0.55%) 1.9506 

105(1/7) 1.8618 (0.46%) 1.8622 (0.48%) 1.8533 

115(3/23) 1.7767 (0.45%) 1.7692 

125(3/25) 1.7020 (0.38%) 1.6956 

135(1/9) 1.6387 (0.50%) 1.6361 (0.34%) 1.6305 



point, while at the same time the JL overlap is only 0.665 
(see Table HTH! 



TABLE V: Values of the ratio a [Eq. O] for JL, REM, 
and exact-diagonalization wave functions for TV = 6 electrons 
and various angular momenta L; v (given in parentheis) are 
the corresponding fractional filling factors. Recall that the 
angular momenta for the JL functions are L 3L — N(N — 
l)(2m + l)/2, m = 0,1,2,3,... 





JL 


REM 


EXACT 


75(1/5) 


5.000 


1.964 


2.877 


105(1/7) 


7.000 


1.972 


2.708 


135(1/9) 


9.000 


1.978 


2.726 



forces^ On the other hand, as we have stressed earlier in 
this paper and in Ref. [lja, the REM functions are able to 
capture the essential effects of the correlations associated 
with the long-range Coulomb force; thus, in agreement 
with the EXD results, the REM a values remain finite as 
L — > oo. 



E. Exponents of current-voltage power law 



SUMMARY 



Another quantity of theoretical and experimental in- 
terest is the ratio 



JL 

max 



1) 



n (t) 



(14) 



of the corresponding occupation numbers at l^ ax — 1 
and The interest in this ratio is due to the fol- 

lowing two facts: (I) The value of a for the JL func- 
tion at different fractional fillings has a particular an- 
alytic value i.e., it is given by a. 3lj (v) = 1/v = 
2m + 1, m = 1, 2, 3, 4, (II) a happens to enter as the 
exponenliSiSi of the voltage in the current-voltage law, 
/ oc V a , for external electron tunneling into an edge of 
a fractional quantum Hall system. Recent investigations 
have found that both the experimental and computed^ 
EXD value of a at v = 1/3 deviates from the JL predic- 
tion of 3, being in all instances somewhat smaller (i.e., 
- 2.7). 

Table |V] displays the values of a for TV = 6 and for 
the JL, REM, and EXD wave functions at various val- 
ues of the total angular momentum L. We have checked 
that our numerical values for ev JL (derived by dividing 
the proper n JL 's; see Fig. [SJ) are equal to 2m + 1 within 
the numerical accuracy. As seen from Table |V] a most 
striking weakness of the JL functions is that the cor- 
responding a JL 's diverge as L — > oo, a behavior which 
contrasts sharply with the EXD values that remain at 
all times finite and somewhat smaller than 3. Such a 
dramatic difference in behavior should be possible to be 
checked experimentally. Furthermore, we note that the 
REM values, although somewhat smaller, they are close 
to the EXD ones and remain bounded as L — > oo. 

We conclude that this dramatic qualitative and quan- 
titative weakness of the JL functions is due to their be- 
ing exact solutions of a family of short range interparticle 



Exact diagonalization (EXD) results for the lowest 
rotational band of a circular QD with TV = 6 elec- 
trons in strong magnetic fields were reported^ here for 
the first time in the range of high angular momenta 
70 < L < 140 (covering the corresponding range of frac- 
tional filling factors 1/5 > v > 1/9). These EXD re- 
sults were used in a thorough assessment of the ability of 
the composite-fermion^/Jastrow-Laughlin 1 - and rotating- 
electron-moleculei^ trial wave functions to approximate 
the exact wave functions in the case of 2D QD's. 

A detailed comparison (addressing five properties, i.e., 
prediction of magic angular momenta, radial electron 
densities, occupation number distributions, overlaps and 
total energies, and exponents of current-voltage power 
law) shows that the REM many-body wave functions 
provide a description that is superior to that obtained 
through the CF/JL ones. An important finding is that 
"global" quantities (like overlaps and total energies) are 
not particularly reliable indices for comparing exact and 
trial wave functions; a reliable decision on the agree- 
ment, or lack of it, between exact and trial wave func- 
tions should include detailed comparisons of quantities 
like radial electron densities and/or occupation number 
distributions. 

We finally note that the CF/JL wave functions have 
been most useful for the modeling of the bulk fractional 
quantum Hall effect. However, the theoretical investiga- 
tions concerning the bulk system have unavoidably, due 
to computational limitations, relied on finite-size systems 
to assess the validity of the CF/JL wave functions. Thus 
it is natural to conjecture that the unexpected finding of 
this paper, i.e., that the CF/JL functions exhibit remark- 
able weaknesses in reproducing the exact wave functions 
of QD's in high B, may have ramifications for our present 
understanding of the fractional quantum Hall effect it- 
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self. Investigations of such probable ramifications, and 
related questions concerning the domain of validity of the 
REM and CF/JL wave functions in the bulk, will be ad- 
dressed in future publications. In the present paper, we 
focused on the case of QD's, which constitute a theoret- 
ically self-contained problem when exact-diagonalization 
calculations become available; in the near future, a wider 



range of such calculations will be within reach, due to 
new generations of powerful computers. 
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tributed uniformly on a parallel disk at a distance d from 
the electron layer (the typical d in experiments is d > 101b, 
see Ref. IHoT) . As was the case with the harmonic external 
potential, these authors found again that their external 
confinement influences which magic-L state becomes the 
ground state of the system. Most importantly, the ground- 
state wave functions in their exact-diagonalization study 
exhibit strong oscillations in the radial electron density [in 
an apparent agreement with the classical ring configura- 
tions of Wigner molecules] and in disagreement with the 
CF/JL wave functions. It is interesting to note the coinci- 
dence, for all practical purposes, of the exact radial elec- 
tron density for TV = 6 and L = 105 calulated by Wan 
et al. with that calculated by us [compare figure 5(d) in 
Ref. |H3 with the middle panel of figure 5 in this paper] . In 
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order to account for the disagreement between the exact 
and CF/JL wave functions, Wan et al. were led to use the 
concept of "edge reconstruction". In the case studied by 
us, however, our exact-diagonalization results (and those 
of Tsiper and Goldman, see Ref. 163) do not include any 
external confinement, a fact that rules out "edge recon- 
struction" as the underlying cause for the disagreement 
between the exact and CF/JL wave functions. As we have 



pointed out in this paper previously (see section IV. E and 
also Ref. IT^l , this disagreement arises from the fact that 
the CF/JL functions do not capture the long-range charac- 
ter of the Coulomb interelectron repulsion. On the contrary 
the REM wave functions are able to capture the long-range 
Coulombic correlations and thus are in better agreement 
with the wave functions from exact diagonalization. 



